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Abstract 

We inspect a particular gauge field theory model that describes the properties of a variety 
of physical systems, including a charge neutral two-component plasma, a Gross-Pitaevskii 
functional of two charged Cooper pair condensates, and a limiting case of the bosonic 
sector in the Salam- Weinberg model. It has been argued that this field theory model 
also admits stable knot-like solitons. Here we produce numerical evidence in support for 
the existence of these solitons, by considering stable axis-symmetric solutions that can 
be thought of as straight twisted vortex lines clamped at the two ends. We compute 
the energy of these solutions as a function of the amount of twist per unit length. The 
result can be described in terms of a energy spectral function. We find that this spectral 
function acquires a minimum which corresponds to a nontrivial twist per unit length, 
strongly suggesting that the model indeed supports stable toroidal solitons. 
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Recently, a gauge field theory model with two charged bosons has been conceived, to describe 
a two-component plasma of negatively and positively charged particles |l| . But this model also 
appears to describe a large variety of other physical phenomena, including a Gross-Pitaevskii 
functional of two band superconductivity and the bosonic sector in the Salam- Weinberg 
model in the limit where the Weinberg angle 6w = |§. In U (see also |§) it has been 
proposed that the model supports stable, self-confining and knot-like solitonic configurations. 
This would be somewhat remarkable, since it would contrast some of the widely held views 
in plasma physics that such configurations of plasma can not exist in general. This is due to 
a simple application of the Shafranov virial theorem which states that a static configuration 
of plasma in isolation is dissipative ||. The proposed model, however, escapes this no-go 
theorem by incorporating non-linear interactions which are not accounted for by mean field 
variables such as the pressure 0. 

The soliton solutions in the gauge model that we shall inspect can be viewed as a bundled 
filaments of twisted magnetic flux lines. The twisting is governed by a certain topological 
quantity, the Hopf invariant. Nontriviality of the Hopf invariant ensures that the flux lines are 
knotted, or linked. Numerical simulations, in the absence of effective analytical tools, seem so 
far to be the best way to help explore the nature of the soliton solutions. But even then their 
intricate topology makes full three dimensional simulations a daunting task. In this letter we 
present and analyse a tractable, yet challenging, simulation of the model, where the magnetic 
flux lines are twisted in an axis-symmetric manner. Such configuration of solutions can then 
be viewed as straight but twisted vortex tubes. 

As such, the stable vortices in the model we study can be applied to study a number of 
interesting physics. They may relate to the coronal loops on the solar photosphere JT|, to 
Meissner effect in two-band superconductors or to higher energy topological configurations 
in the weak sector of the standard model 0, ||. Our study will then serve as a test bed to 
understand knot solitons in general. Serious, three dimensional searches for knotted structures 
in field theory models are attempted only recently. One prototype model, initiating these 
searches 0, [§ (see also 0, [10], ITTJ), is a Skyrme like 0(3) non-linear sigma model. This 
model could be envisaged as describing the infrared phase of the pure SU(2) Yang-Mills 
theory, with glueballs represented by the knotted flux tubes of gluon RI2| . 
We start from the classical kinetic theory model of a two-component plasma of electromag- 
netically interacting electrons and ions, given by the non-relativistic action [|TJ, 
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Here, \l/ e and are the two complex non-relativistic, macroscopic Hartree wave functions 
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describing the electrons (e) and ions 
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i) with their respective masses m and M. Numerically, 
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The electron and ion densities are, respectively, 
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given by \l/*\l/ e and and their total integrals over the three-space give the total electron 

number N e and the total ion number iVj. Charge neutrality requires N e = iVj. 
The ensuing static Hamiltonian in the Coulomb gauge is, 
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where B is the magnetic field. We note the similarity of the above with the Hamiltonian that 
describes the the bosonic sector of the Salam- Weinberg model at 9w = 0: At this prescribed 
value of the Weinberg angle the masses of and Z boson become infinite and they, therefore, 
decouple from the theory. Now, assigning the hypercharge matrix of the Higgs doublet to be 
proportional to the third Pauli matrix, t$, the static Hamiltonian of the bosonic sector of the 
Salam- Weinberg model is, 
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where the Higgs doublet is given by, 




In the limit of weak self-couplings between the Higgs fields H$w is notably similar to the 
Hamiltonian in Eqn.(Q) with the obvious identification <p +, ~ = 

An effective static energy functional of plasma can be obtained from Eqn.(^j) in a self-consistent 
gradient expansion by keeping terms with at most fourth order in the derivatives of the vari- 
ables. In order to describe the ensuing tubular field configurations appearing in the model, it 
is natural to introduce a new set of variables ill, 
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where a is a parameter expressed in terms of the reduced mass, /i, through the relation 
fi = msin 2 a = Mcos 2 a, and p 2 is related to the plasma density, and 9 is a shape function 
measuring roughly the distance from the center of the configuration, and ip and x are the 
toroidal and poloidal coordinates in R 3 . By defining a three component unit vector, 



if = (cos (x + <p) sin 9, sin (x + <p) sin 9, cos i 
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it can be shown that the static energy is ||, 

E = j d^lC^dkpf + C 2 p 2 \d k H\ 2 + c z {v: -drfxd^f (5) 
+ C 4 p 4 (cos 2a — cos#) 2 ] , 
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where C\ = — sin 2 2a, C 2 = — = sin 2 2a, C 3 = — — , and C 4 = -. The effective cou- 
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pling g describes the remnant of the Coulomb interaction in the plasma, in the limit of short 
Debye screening length. At this point, it is of interest to compare the above energy density 
with that of ||. There, the energy density consists of the two middle terms in the above 
expression, for a constant p. Indeed, the presence of a nontrivial coupling between p and if 
in the above expression for energy density is especially noticeable. 

In order to have finite energy configurations, asymptotically at large distances if must go to a 
constant value with n 3 = cos 2a, and also p = Po asymptotically at large distances. Here, po is 
a constant valued characteristic plasma parameter related to the plasma density at the bulk. 
For example on the solar photosphere po is of order of magnitude 10 15 /m 3 . The unit vector 
if, when combined with the boundary conditions, describes a map from the compactified R 3 
to the target S 2 . Under this map the pre- image of a point on the target is generically a circle, 
knotted or linked, and such circle is a constituent element of the magnetic field lines in the 
plasma. Any two pre-image circles are linked with their Gauss linking number given by the 
topologically invariant, integer valued Hopf number H, 

Stable finite energy knotted and linked soliton solutions are classified by the Hopf invariant 
H . A non-trivial question of interest is to answer as to for which Hopf numbers the solutions 
are actually knotted, not merely linked. 

The equations of motion arising from varying Eqn.([5D depend on two parameters po and 
g. However, by re-scaling p — ► pop and x — > xqx, where p and x are both dimensionless 
quantities, the equations of motion can be recast to make dependent only on g. Henceforth, 
all the expressions are written in terms of the dimensionless variables p and x and we continue 
to denote them as p and x, respectively. The parameter xq has the dimension of length and 
has the expression xl = 3 2 . 

We are interested in the axially symmetric solutions to obtain straight twisted vortices. The 
vector field generating such axis-symmetric twist is 

V = (ify - \d z ) 

and the Lie derivative of the field variables with respect to V must be zero. The ansatz for 
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the fields, satisfying the previous condition, are: (x + f) = kz + ar, p = p (r), and 9 = 9 (r). 
Here, k is a real number, and a denotes the twist per unit length. Similar ansatz were used 



in pT5fl . We consider the tube to be clamped at its two ends and the length of the tube is L. 

kaL 

The Hopf invariant now becomes H = — — . Notice that, had we taken the tube length to 
be infinite, the Hopf number would have become infinite. It might be tempting to consider 
a straight tube with the topology of a torus, but this does not work: The toroidal topology 
implies that the fields are periodic in z with period L, which in turn implies that kL is an 
integer multiple of 2tt. One would then have a to be a rational number only. This is why we 
exclude henceforth a straight tube with only toroidal topology. 
The energy functional Eqn.(BI) in the axially symmetric ansatz reads, 
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and the Coulomb coefficient C 
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of the parameter dependence of the field variables of the above energy functional is particularly 
revealing, as we will see shortly. 

The numerical solution to the Euler equations of motion arising from (^) are obtained by 
seeking the fixed points of the following system of equations: 
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The simulations are run on a lattice of finite size. At one boundary end, we take p = 1, and 
9 = 2a. At the other end, which is the origin, we also need the boundary values of p and 9. 
For these values to be fixed, first note that the equations of motions are invariant under the 
parity transformation r — > — r and, therefore, we are free to choose independently 9 and p to 
be either odd or even function in r. We have 9(0) = corresponding to choosing 9 to be odd. 
As we want p to be non-vanishing at the origin, we take p to be even. It should be mentioned 
that p could have been chosen equally well to be an odd function, but we would like charge 
densities to be non-vanishing at the origin. With the boundary conditions so chosen and for 
technical reasons in order to facilitate the simulations, the range of the lattice is also extended 
to the negative values of r. Next, we choose initial profiles for p and 9 matching with the 
boundary conditions. Finally, the equations are solved for fixed k — 1, since higher k would 
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corresponds to the configurations with higher energies and are, therefore, excluded from our 
simulations. In the simulations we have performed the Coulomb term C is chosen to be of 
small value, 0.1, 1.0, and 10.0, and the twist parameter a is made to lie in the range [0.2, 2.0]. 
In Figure 1. we have drawn plots for the energy per unit Hopf number, S/H, against the 
twist per unit length a. Each point on the plot corresponds to a solution of the equations of 
motion for a given a. These energy plots for different Coulomb couplings can be described 
by spectral functions f(a,C). As visible from the plots, these spectral functions all have the 
following features in common: For each C the spectral functions are smooth, positive, and 
strictly convex with a nontrivial minimum. That for a given C the spectral function f(a, C) 
is a positive convex function of the twist a could be seen apriori from the form of the energy 
Eqn.(0). As both for a — > and oo, S/H diverges, given that the solutions are smooth. 
However, what is remarkable is the form of the graph of the function. The unique minimum 
point of the graph, occurring at ac, represents the true stable solution that an axis-symmetric 
vortex tube with a given number of twist would settle to. Too many, or too few, twists per 
unit length in the tube to begin with would result in instability. It is to note that as C 
varies so does ac, but little. We have performed the numerical simulations for a number of 
representative values of the parameters, quite far away from realistic values applicable e.g. to 
coronal loops on the solar photosphere. It would certainly be of interest to extrapolate our 
calculations for the physically interesting values of C . Unfortunately, in this case the various 
numerical parameters involved deviate from each other by several orders of magnitude. As 
a consequence we find numerical intractability as a hindrance for achieving this goal, and 
postpone it to future publications. 

The plots for the energy and ion densities have also some interesting features, as described in 
Figs. 2 and 3. Namely the peak of the energy density plot lies, somewhat counter-intuitively, 
slightly off the center. The reason for this can be traced to the twisting of the field lines. By 
looking at the energy density plot, one can furhtermore estimate the thickness of the vortex 
tube and, on the other hand, by reading off the minimum point of the spectral curve one can 
estimate its length. We find that the ratio between the length and the thickness turns out to 
be 2.5. This suggest that for a would-be toroidal configuration the energy density becomes 
lumped at the center of the toroidal structure in analogy with the model |7j ; see [|1(], [T]J . 
The total number of ion and electron numbers, respectively, N e and iVj, are tabulated in Table 
1 for different values of a and C. Clearly, we do not obtain N e = iVj implying the vortices 
carry charge, but this we consider to be a finite size-effect as the simulations are run on a finite 
lattice. One can, however, conclude from the table that the heavier ions are concentrated more 
towards the center of the tube and the lighter electrons are spread out more to the bulk. 
To conclude, we have presented numerical evidence that the gauge field theory model of plasma 
|U does admit stable knotted or linked soliton solutions. We have searched for a particular 
kind of soliton in the shape of straightened twisted tube. The length of the tube is clamped 
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at a fixed length and the stability of the solution requires the tube to be twisted by a certain 
amount per unit length. It would be of interest to extend our analysis to the toroidal case, 
where the addition of curvature demands a full three dimensional simulation. And it would 
certainly be desirable to run simulations in the physically interesting regime of the param- 
eter space in order to understand for example the origin of the coronal loops on the solar 
photosphere, or knotted configurations in the Salam- Weinberg model. At the moment this 
is hampered by technical reasons as the various parametes involved deviate from each other 
by several orders of magnitude, in particular in the case of solar surface. This unfortunately 
undermines the numerical stability in our present simulations. Besides finding plausible appli- 
cations in areas of plasma and condensed matter physics, our study appeals also directly to the 
understanding of knot solitons in gauge field theories in general. Indeed, a highly interesting 
question would be whether the weak sector of the standard model admits knot solitons ||. |6j . 
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Table 1: Total number of ions, Ni, and electrons, N e , for different values of the twist per unit 
length, a, and the Coulomb coupling, C. 
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Figure 1: The total Energy per unit Hopf number as a function of a, for different values of C. 
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Figure 2: The energy density versus distance for different values of C. 
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Figure 3: Ion density versus distance at the minimum ac for different values of C. 
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